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ABSTRACT 


The  problem  of  pointing  an  instrument  on  the  earth's  surface  at  an 
object  in  orbit  about  the  earth  is  investigated.  Necessary  expressions  are 
derived  for  calculating  continuously  azimuth,  elevation,  range,  and  range 
rate  from  inputs  of  instrument  site  and  orbit  parameters.  An  electro¬ 
mechanical  orbit  computer  capable  of  performing  the  required  mathematical 
operations  is  described  Test  results  are  listed  for  a  similar  computer 
used  against  a  typical  orbiting  body. 
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GLOSSARY 


A  elevation  angle.  The  angle  between  the  local  horizontal  plane  and 

the  line  of  sight  measured  in  the  vertical  plane  through  the  line  of 
sight. 

a  semi-major  axis.  One-half  the  length  of  the  major  axis  of  the 

ellipse  of  orbit 

B  azimuth  angle.  The  dihedral  angle  between  the  plane  of  the  north 

meridian  at  the  site  and  the  vertical  plane  through  the  line  of 
sight,  measured  about  the  site  zenith  in  a  clockwise  direction 
from  north, 

b  semi-minor  axis  One-half  the  length  of  the  minor  axis  of  the 

ellipse  of  orbit. 

c  focal  distance.  The  distance  between  the  center  of  the  ellipse  of 

orbit  and  a  focus  of  the  ellipse. 

E  eccentric  anomaly.  The  angle  at  the  center  of  the  ellipse  between 

the  radius  to  the  point  of  perigee  and  the  radius  to  the  orthogonal 
projection,  from  the  major  axis,  of  the  target  on  the  circum¬ 
scribed  circle,  measured  from  perigee  in  the  direction  of  orbital 
motion 

e  eccentricity  The  eccentricity  of  the  ellipse  of  orbit,  defined  as 

the  ratio  between  the  focal  distance  and  the  length  of  the  semi¬ 
major  axis,  or  c/a. 

f  flatness.  The  flatness  of  the  earth's  surface,  defined  as  the  ratio 

between  the  difference  in  semi-major  and  semi-minor  axis  lengths 
and  the  length  of  the  semi-major  axis.  Its  presently  accepted 
value  is  1/297.0. 

h  slant  range  The  distance  from  the  site  to  the  target,  measured 

along  the  line  of  sight. 

inclination.  The  dihedral  angle  between  the  plane  of  the  celestial 
equator  and  the  orbital  plane.  It  is  measured  about  the  ascending 
node  counterclockwise  from  the  equator  to  the  first  intersection 
with  the  orbit  plane. 


1 


LST 


local  sidereal  time.  The  dihedral  angle  between  the  plane  of  the 
site  meridian  and  the  plane  of  the  meridian  through  the  first  point 
of  Aries.  It  is  measured  about  the  north  pole  from  the  first  point 
of  Aries  toward  the  east 

M  mean  anomaly.  The  angle  at  the  center  of  the  earth  between  the 

radius  to  the  point  of  perigee  and  the  radius  to  a  point  the  target 
would  occupy  at  the  time  under  consideration  if  it  were  in  a 
circular  orbit  about  the  earth.  This  angle  is  equal  to  the  time 
since  perigee  passage  multiplied  by  the  average  orbital  rate,  or 
M  =  2tt/P  (LST  -  T).  It  is  measured  from  perigee  in  the  direction 
of  target  motion, 

P  period  of  orbit.  The  time  required  for  the  target  to  complete  one 

orbit. 

R  radius  of  the  earth.  The  distance  from  the  center  of  the  earth  to 

the  site. 

r  geocentric  range.  The  distance  from  the  center  of  the  earth  to  the 

target. 

S  orbit  area.  The  area  of  the  ellipse  of  orbit. 

V  range  angle.  The  angle  at  the  center  of  the  earth  between  the 

radius  to  the  ascending  node  and  the  radius  to  the  target,  measured 
from  the  ascending  node  in  the  direction  of  target  motion  along  the 
orbit. 

range  angle  of  perigee.  The  angle  at  the  center  of  the  earth 
between  the  radius  to  the  ascending  node  and  the  radius  to  the 
point  of  perigee,  measured  from  the  ascending  node  in  the 
direction  of  target  motion  along  the  orbit. 

5  latitude  correction.  The  angular  difference  between  geodetic 

latitude  X.  and  geocentric  latitude  X  . 

g  ° 

X  geodetic  latitude.  The  angle  between  the  equatorial  plane  and  the 

perpendicular  to  the  earth's  surface  at  the  site,  measured  in  the 
plane  of  the  meridian  through  the  site. 

geocentric  latitude.  The  angle  between  the  equatorial  plane  and 
the  radius  from  the  center  of  the  earth  to  the  site,  measured  in 
the  plane  of  the  site  meridian. 
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perigee  passage.  The  time  at  which  the  target  passes  perigee. 


Q,  longitude  of  the  ascending  node.  The  dihedral  angle  between  the 

plane  of  the  ascending  node  meridian  and  the  plane  of  the  first 
point  of  Aries  meridian,  measured  about  the  north  pole  in  the  east 
direction  from  the  first  point  of  Aries. 
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I. 


INTRODUCTION 


Location  of  a  satellite  in  unperturbed  orbit  about  the  earth  may  be 
described  in  terms  of  six  parameters  Location  of  a  site  on  the  surface  of 


the  earth  may  be  described  in  terms  of  three  parameters.  From  these 
parameters  continuous  pointing  angles,  range,  and  range  rate  may  be 
calculated.  A  special  effort  is  made  in  this  report  to  present  the  analytic 
solution  to  the  orbital  pointing  problem  in  such  a  manner  that  the  suggested 
machine  solution  is  self-evident. 


II 


ANALYTIC  SOLUTION 


When  target  position  is  known  in  its  plane  of  orbit,  it  is  possible  by 


a  series  of  successive  rotations  and  translations,  determined  by  the  orbit 
and  site  parameters,  to  arrive  at  the  coordinates  of  target  position  with 
respect  to  a  site  of  known  location.  Assume  from  Fig.  1  that  the  range 


angle  from  perigee  (  y  -  y  )  is  known  Then 

x  =  r  sin  (y  -  y  ) 

1  o 


y =  r  cos  ( y  -  y  ) 


1  o 


z 


1 


0 


(1) 


The  coordinate  system  is  rotated  about  its  Z^-axis  through  the  angle  y  so 
that  the  Y^-axis  passes  through  the  ascending  node.  In  the  new  system 


x  =  x  cos  y  f  y  sin  y 
Z  1  o  1  o 


y_  =  y.  cos  v  -  x  sin  V 


2  3  1  o  1  o 


(2) 


Next,  the  system  is  rotated  about  the  Y^-axis  through  the  inclination 
angle  i,  as  shown  in  Fig.  Z.  The  X^-axis  is  in  the  equatorial  plane  and 
passes  through  the  meridian  90°  east  of  the  ascending  node;  the  Y^-axis  is 
in  the  equatorial  plane  and  passes  through  the  ascending  node;  the  Z^-axis 


-axis 
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passes  through  the  north  pole.  Target  coordinates  are 

cos  i  , 

y3  =  y2  . 

Z3  =  X2  S^n  *  ^ 

Then,  a  rotation  about  the  Z^-axis  through  the  angle  (FST  -ft')  places  the 
Y^-axis  through  the  meridian  of  the  site 

x^  =  x^  cos  (LST  -ft)  -  y  sin  (LST  -  £2)  , 

cos  (LST  -  ft)  +  sin  (LST  -  ft)  , 

z  =  z  .  (4) 

4  3 


Next,  the  coordinate  system  is  rotated  about  the  X^-axis  through  the  angle 

(90  -  X  )  to  place  the  Z  -axis  perpendicular  to  the  horizontal  plane  at  the 

g  5 

site,  as  shown  in  Fig  3. 

xc  =  x 
5  4 


y^  =  y  .  sin  ^  -  z  ,  cos 

5  4  g  4 


X 


z  =  z  sin  X  +  y  cos  X 
5  4  g  4 


(5) 


Translations  along  the  Y  -  and  Z  -axes  are  required  to  position  the 

5  5 

origin  of  the  coordinate  system  at  the  antenna  site.  From  Fig.  3  the 

angular  difference  6  between  geodetic  latitude  X  and  geocentric  latitude  X 

g  ^ 

is  calculated  The  slope  of  the  perpendicular  to  the  surface  at  (y^,  z^)  1S 


(1  -  f)‘ 


where  f  is  the  flatness  of  the  earth's  surface.  The  slope  of  the  radius  to 
the  same  point  is  z^/y^  The  difference  between  the  two  latitudes 


2 


„  {2f  -  f  )  sm  X  cos  X 

6  =  tan"  - 2 - *2 - 2^ 

cos  X  h  ( 1  -  f )  sin  X 


or,  since  6  is  small  and  f  is  much  less  than  unity, 


6  ~2f  sin  X  cos  X  -  f  sm  2  X 
g  g  g 


The  origin  of  the  coordinate  system,  then,  is  translated  to  the  site  by 
making 


X6 

=  X5 

> 

y6 

=  y5 

-  R  sin 

5  s 

ay  -  Rf  s in  2  X. 
5 

z  . 

-  z 

-  R  cos 

6  i 

-  R 

6 

5 

5 

The  X  -axis  is  horizontal  in  the  east  direction,  the  Y  -axis  is  horizontal 
o  o 

in  the  south  direction,  and  the  Z^-axis  is  vertical.  From  Fig.  4  the 
pointing  angles  are  for  azimuth 


B  -  tan 


-1 


-Yi 


and  for  elevation 


A  =  tan 


-1 


,  2  21/2 
+  (x6  f  y6) 


(7) 


(8) 


The  slant  range  is 

222  1/2 

h  =  (x6  f  y6  +  z6)  .  (9) 

Range  rate  is  the  time  derivative  of  Eq.  (9). 

To  calculate  pointing  information  in  real  time,  it  is  necessary  to 
establish  a  relationship  between  time  and  y  -  y  ,  the  orbital  range  angle 
measured  from  perigee  to  the  target.  Referring  again  to  Fig.  1,  other 
expressions  for  the  components  of  target  position  are 
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(10) 


=  b  sin  E  , 

Y1  =  a<cos  E  -  e) 

From  these  expressions,  range  from  the  center  of  the  earth  to  the  target 
is 


r  -  (x  +  y  )  =  a(l  -  e  cos  E) 

Settings  Eqs.  (10)  equal  to  the  corresponding  Eqs.  (1) 

r  cos  (v  -  Y  )  =  a(cos  E  -  e) 
o 


(11) 


2  1/2 

r  sin  (\  -  V  )  =  b  sin  E  =  a(l  -  e  )  sin  E 
o 


(12) 


F rom  Eqs .  (12) 


d(v  -  v  ) 
o 


M  V/2 

(1  -  e  ) 

1  -  e  cos  E 


dE 


(13) 


From  Fig.  1  the  incremental  area 


2  2  2  1/2 
dS  =  A—  d(Y  -  Y  )  =  At-  (1  -  e  )  '  (1  -  e  cos  E)dE 
2  o  2 


(14) 


Kepler’s  second  law  for  orbiting  bodies  states  that  the  time  derivative  of 
the  area  swept  out  by  the  radius  (r)  from  a  focus  of  the  ellipse  to  the 
orbiting  body  is  a  constant.  Making  use  of  this  principle, 


dS 

dt 


=  K  = 


area 

period 


Trab 


Tra2(l 


2.1/2 

e  ) 


(15) 


and 


dS 


Tra2(l 


2.1/2 

e  ) 


dt 


Combining  Eqs.  (14)  and  (16), 


2tt 

P~ 


dt  =  ( 1  -  e  cos  E)dE 


(16) 


(17) 
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Integrating  both  sides  of  Eq.  (17)  over  the  time  elapsed  since  passage  of 
perigee , 


( 1  -  e  cos  E)dE 


(18) 


2tt 


(LST  -  T)  =  E  -  e  sin  E 


P 


or 


M  =  E  -  e  sin  E 


(19) 


From  Eq.  (19)  the  eccentric  anomaly  E  may  be  determined  as  a  function  of 
time.  The  range  angle  from  perigee  ( y  -  y  )  is  calculated  by  Eqs.  (12) 
after  E  is  determined. 


III.  MACHINE  SOLUTION 


The  series  of  operations  necessary  to  proceed  from  inputs  of  site 


and  orbit  parameters  to  the  desired  outputs  is  accomplished  readily  by  an 
electromechanical  analog  computer.  A  functional  diagram  of  such  a 
machine  is  shown  in  Eig.  5. 

The  basic  component  used  in  the  computer  is  the  resolver,  which  is 
capable  of  performing  almost  all  the  mathematical  manipulations  required, 
It  accepts  two  electrical  inputs,  u^  and  v^  ,  and  a  mechanical  angle  input 
Its  two  electrical  outputs  are 


(20) 


v  =  v  cos  \\>  +  u  sin  1 

L  1 


which  describe  the  operations  necessary  for  rotating  a  rectangular 
coordinate  system.  Such  an  operation  is  required  for  Eqs.  (1),  (2),  (3), 


(4),  and  (5).  If  components  u^  and  v  are  given  and  the  angle  \\>  is  driven 
until  u^  in  Eq,  (20)  is  equal  to  zero,  then 


=  tan 


5 


(21) 


,2  21/2 

V2  -  (U1  +  V 


This  technique  is  used  to  solve  Eqs.  (7),  (8),  and  (9).  The  only  other 
necessary  operation  is  algebraic  summing,  which  is  performed  by  parallel 
input  operational  amplifiers. 

From  Fig.  5,  one  of  the  required  inputs  to  the  computer  is  AM, 

the  time  integral  of  the  mean  anomaly  rate.  This  shaft  rotation  is  produced 

by  the  mechanism  shown  in  Fig  6.  A  synchronous  motor  is  driven  at 

constant  speed  by  a  constant  frequency  source.  The  motor  drives  through 

reduction  gearing  the  inputs  to  23  clutches  which  are  capable  of  switching 

23  different  speeds  into  the  output.  The  output  rate  is  the  sum  of  the  output 

rates  of  the  clutches  which  are  engaged  Gearing  is  so  arranged  that  the 

relative  speeds  of  the  clutch  outputs,  starting  at  the  least  significant  clutch, 

are  1,  2,  4,  8,  10,  20,  40,  80,  etc.  Thus,  five  full  decades  and  one  partial 

decade  of  discreet  rates  may  be  selected  by  controlling  the  clutches  in 

-  5 

8-4-2- 1  binary  coded  decimal  format.  Rates  from  1  x  10  degree  per 

-5 

minute  to  8.  0  degrees  per  minute  are  available  in  increments  of  1  x  10 
degree  per  minute.  Accuracy  of  the  average  output  rate  is  dependent  only 
upon  the  accuracy  of  the  constant  frequency  source. 

The  computer  requires  a  total  of  nine  inputs. 


1. 

e 

eccentricity  of  the  orbit. 

2. 

R/a 

ratio  between  radius  of  earth  at  site  and 

semi-major  axis  of  orbit. 

3. 

M 

o 

initial  value  of  mean  anomaly. 

4. 

dM/dt 

mean  anomaly  rate. 

5. 

range  angle  of  perigee 

6. 

i 

orbit  inclination. 

7. 

(LST) 

o 

initial  value  of  local  sidereal  time. 

8. 

Q 

longitude  of  ascending  node. 

9. 

g 

geodetic  latitude  of  site. 
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Changes  to  mean  anomaly  M  and  local  sidereal  time  LST  are  generated  by 
the  machine  to  provide  continuous  solutions  in  real  time. 

IV.  IMPLEMENTATION 

Accuracy  of  the  computer  depends  to  a  great  extent  upon  the  accuracy 
of  the  resolvers  used  The  selected  winding  compensated  resolvers  have  a 
functional  error  of  less  than  0  05  percent  of  output  at  maximum  coupling 
and  total  interaxis  error  of  less  than  1  5  minutes  of  arc.  Gain  of  unity  and 
phase  shift  of  zero  between  resolver  inputs  and  resolver  outputs  is  main¬ 
tained  by  making  use  of  high  gain  solid  state  buffer  amplifiers.  The 
resolver  compensator  winding  is  magnetically  coupled  to  the  primary 
winding  so  that  it  sees  any  changes  in  transformation  ratio  between  the 
primary  and  secondary  windings.  The  input  voltage  and  compensator  wind¬ 
ing  voltage  are  compared  at  the  input  of  the  buffer  amplifier  Any  change 
in  the  resolver  transformation  characteristics  will  be  reflected  by  a  change 
in  compensator  voltage.  This,  in  turn,  will  adjust  the  closed  loop  gain  of 
the  resolver-buffer  amplifier  combination  to  restore  transformation  ratio 
to  its  original  value  The  high  input  impedance  of  the  buffer  amplifier  is 
also  useful  in  providing  isolation  between  interconnected  resolvers. 

A  number  of  parameters  must  be  introduced  into  the  computer  by 
accurate  angle  setting  of  resolver  shafts  A  separate  dial  system  and 
associated  precision  reduction  gearing  for  each  resolver  is  costly  and  space 
comsuming.  If,  instead,  a  dual  speed  multipole  synchro  is  coupled  directly 
to  the  resolver  shaft,  it  will  provide  indication  of  shaft  angle  to  within  20 
seconds  of  arc.  A  single  precision  dial  system  with  another  multipole 
synchro  is  used  to  set  consecutively  each  resolver  shaft  by  servo  slaving. 

The  complete  orbit  computer  package  will  occupy  about  16  cubic  feet 
and  weigh  about  500  pounds. 

V  TEST  RESULTS 

Pointing  data  from  a  computer  similar  to  that  described  have  been 
compared  to  the  calculated  positions  of  a  satellite  in  a  medium  altitude 
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orbit  about  the  earth  Observed  errors  for  a  representative  group  of 
satellite  passes  are  tabulated  below. 


^Total  Pointing 
Error 
(Degrees) 

Range 

Error 

(Nautical  Miles) 

Range  Rate 

Error 

(Nautical  Miles/Minute) 

Average 

0.  091 

9.  2 

1.  4 

RMS 

0.  100 

11.  0 

1.  9 

Maximum 

0.  230 

25.  0 

3.  5 

The  magnitudes  of  these  errors  are  in  close  agreement  with  those  calculated 
from  the  known  tolerances  of  the  various  components  used  in  the  computer, 

VI.  CONCLUSIONS 

The  computer  described  is  capable  of  the  simple,  straight-forward, 
and  economical  solution  of  the  orbital  pointing  problem.  It  may  be  used, 
as  well,  against  targets  which  appear  fixed  in  space  and  against  orbiting 
belts.  Well-known  techniques  and  conventional,  readily  available  compo¬ 
nents  are  used  Solutions  are  instantaneous,  continuous,  and  in  real  time 
The  accuracy  obtained  appears  to  be  compatible  to  that  with  which  the  orbit 
parameters  of  most  earth  satellites  are  known. 


*Total  pointing  error  is  the  great  circle  distance  between  the  target  and 
the  line  of  sight.  It  is  the  resultant  error  due  to  the  individual  errors 
about  the  two  pointing  axes. 
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Fig.  2.  Location  of  target  and  site  in  coordinate  system. 
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3-76-1884 
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3-76-1885 


Fig.  4.  Location  of  target  in  site -centered  coordinate  system. 
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rig.  5*  Functional  diagram  -  analog  orbit  computer 
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Fig,  6,  Mean  anomaly  and  local  sidereal  time  driver 


